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Abstract 

The timelike geodesic equations resulting from the Kerr gravitational 
metric element are derived and solved exactly including the contribution 
from the cosmological constant. The geodesic equations are derived, by 
solving the Hamilton- Jacobi partial differential equation by separation of 
variables. The solutions can be applied in the investigation of the motion 
of a test particle in the Kerr and Kerr-(anti) de Sitter gravitational fields. 
In particular, we apply the exact solutions of the timelike geodesies i) to 
the precise calculation of dragging (Lense-Thirring effect) of a satellite's 
spherical polar orbit in the gravitational field of Earth assuming Kerr 
geometry, ii) assuming the galactic centre is a rotating black hole we 
calculate the precise dragging of a stellar polar orbit aroung the galactic 
centre for various values of the Kerr parameter including those supported 
by recent observations. The exact solution of non-spherical geodesies in 
Kerr geometry is obtained by using the transformation theory of elliptic 
functions. The exact solution of spherical polar geodesies with a nonzero 
cosmological constant can be expressed in terms of Abelian modular theta 
functions that solve the corresponding Jacobi's inversion problem. 
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1 Introduction 



1.1 Motivation 



Most of the celestial bodies deviate very little from spherical symmetry, and 
the Schwarzschild spacetime is an appropriate approximation for their gravita- 
tional field [2]. However, for some astrophysical bodies the rotation of the mass 
distribution cannot be neglected. A more general spacetime solution of the grav- 
itational field equations should take this property into account. In this respect, 
the Kerr solution [7] represents, the curved spacetime geometry surrounding a 
rotating mass [8]. Moreover, the above solution is also important for probing 
the strong field regime of general relativity [9]. This is significant, since general 
relativity has triumphed in large scale cosmology [5, 4, 51, 10], and in predicting 
solar system effects like the perihelion precession of Mercury with a very high 
precision [1, 3]. In a recent paper [3] Kraniotis and Whitehouse, provided the 
exact solution in closed analytic form, of the time-like geodesies that describe 
the motion of a test particle in the Schwarzschild gravitational field including 
the contribution from the cosmological constant. The solution for the orbit was 
expressed in terms of Abelian hyperelliptic genus 2 modular functions using the 
solution of Jacobi's inversion problem for hyperelliptic integrals. For zero cos- 
mological constant the solution was provided by the Weierstrafi Jacobi modular 
form. Subsequently, the authors applied the solution to the compact calculation 
of the orbit and the perihelion precession of Mercury around the Sun and com- 
pared with current experimental data. The impressive agreement of the precise 
theory developed in [3] for the perihelion precession and the eccentricity of the 
orbit, based on the exact solution of the geodesic equations, with experiment, 
the precise determination of the cosmological constant effect in conjuction with 
the dedicated efforts of the experimentalists to measure the orbit of Mercury 
with higher accuracy in a series of experiments (BcpiColombo ES A [37] , Mes- 
senger NASA [38]) motivates the developement of the theory beyond the static 
nature of the Schwarzschild spacetime. In this way, one can determine all the 
relativistic effects, and compare the theory with further measurements of the 
cosmological and orbital physical parameters. 

Furthermore, as was discussed in [3], the investigation of spacetime struc- 
tures near strong gravitational sources, like neutron stars or candidate black 
hole (BH) systems is of paramount importance for testing the predictions of the 
theory in the strong field regime. The study of geodesies are crucial in this re- 
spect, in providing information of the structure of spacetime in the strong-field 
limit. The generalization of the precise theory, to the gravitational field of a 
rotating mass, will have important applications also in this domain. 

The study of the geodesies from the Kerr metric are additionally motivated 
by recent observational evidence of stellar orbits around the galactic centre, 
which indicates that the spacetime surrounding the Sgr A* radio source, which is 
believed to be a supermassive black hole of 3.6 million solar masses, is described 
by the Kerr solution rather than the Schwarzschild solution, with the Kerr 
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parameter [22] 

n J -, = 0-52 (±0.1, ±0.08, ±0.08) (1) 

where the reported high-resolution infrared observations of Sgr A* revealed 'qui- 
escent' emission and several flares. This is half the maximum value for a Kerr 
black hole [11]. In the above equation J 1 denotes the angular momentum of 
the black hole (The error estimates here the uncertainties in the period, black 
hole mass (Mbh) and distance to the galactic centre, respectively; G is the 
gravitational constant and c the velocity of light.) 

Taking into account the cosmological constant A contribution, the general- 
ization of the Kerr solution is described by the Kerr -de Sitter metric element 
which in Boycr-Lindquist (BL) coordinates 2 is given by [15, 16]: 



where 



A 2 2 

ds 2 = (cdt - a sin 2 8d<f>) 2 - ?—dr 2 - ^-d9 2 

E 2 p 2 v YJ A r Ag 

A e sin 2 



(aedt - (r 2 + a 2 )d(t>) (2) 



A 2w 2 2^ IGMr 

(l--r 2 )(r 2 +a 2 ) 



3 /v ' 6 



.2 



a 2 A 2 n 
Ag := 1 + — cos 2 6 



1 + P 2 :=r 2 + a 2 cos 2 l 



(3) 



The above solution is stationary, but not static. 

It is the purpose of this paper, to derive the geodesic equations that describe 
the motion for a test particle in the Kerr and Kerr -de Sitter gravitational fields, 
and obtain exact solutions of the corresponding equations in an interesting class 
of possible types of motion, generalizing the results of [3] . 

We also apply the exact solutions of the geodesic equations to the following 
situations: 

Frame dragging from rotating gravitational mass: An essential property of 
the geodesies in Schwarzschild spacetime is that although the orbit precesses 
relativistically it remains in the same plane; the Kerr rotation adds longitudinal 
dragging to this precession. For instance, in the spherical polar orbits we will dis- 
cuss in the main text, (where the particle traverses all latitudes, passes through 

1 J = ca where a is the Kerr parameter. The interpretation of ca as the angular momentum 
per unit mass was first given by Boyer and Price [13]. In fact, by comparing with the Lense- 
Thirring calculations [14] they determined the Kerr parameter to be: a = — 2 ™ , where S7 
and I denote the angular velocity and radius of the rotating sphere. 

2 These coordinates have the advantage that reduce to the Schwarzschild solution with a 
cosmological constant in the limit a — > 0, see [12] . 
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the symmetry axis z, infinitely many times) the angle of longitude increases 
after a complete oscillation in latitude. This phenomenon, is in accordance with 
Mach's principle. 

We shall calculate the dragging of inertial frames in the following situations, 
a) Dragging of a satellite's spherical polar orbit in the gravitational field of 
Earth assuming Kerr geometry, b) Dragging of a stellar, spherical polar orbit, 
in the gravitational field of a rotating galactic black hole. 

The material of this paper is organized as follows. In section 2, starting with 
the Kerr metric we derive the geodesic equations by integrating the Hamilton- 
Jacobi partial differential equation by separation of variables. In section 3, we 
solve exactly the timelike spherical polar geodesies. The solution for the orbit is 
expressed in terms of the Weierstrafi elliptic function. We then apply the exact 
solution in detcrming the frame-dragging (Lcnsc-Thirring effect) of a satellite's 
spherical polar orbit in the gravitational field of Earth as well as the dragging of 
a stellar spherical polar orbit by a galactic black hole. The exact expression of 
dragging is proportional to the real half-period of the Weierstrafi Jacobi mod- 
ular form. An alternative expression in terms of a hypergeometric function is 
also provided. In section 4, we solve exactly spherical non-polar geodesies. In 
section 5, we describe the general solution of (non-spherical) timclike geodesies 
in Kerr metric, using the transformation theory of elliptic functions, an ap- 
proach developed by Abel in [17]. In section 6, we derive the timelike geodesic 
equations in Kerr metric that includes the contribution from the cosmological 
constant. The derivation is produced, by integrating the Hamilton- Jacobi par- 
tial differential equation again by a separation of variables. We subsequently 
discuss the exact solution, for spherical polar geodesies with a cosmological con- 
stant. The integration can be achieved, by solving the corresponding Jacobi's 
inversion problem for a system of Abclian integrals that we have formulated. 
Afterwards, the equatorial geodesis with a cosmological constant are presented 
together with the exact solution of the corresponding circular equatorial orbits. 
The application of the exact solutions with the A term present will be a subject 
of a future publication. Finally, in section 7 we present our conclusions and 
discuss future research directions. In the appendices, we collect some of our 
formal calculations as well as some useful formulas for the solution of cubic and 
quartic equations and definitions of genus-2 theta functions. 

2 Geodesies and the separability of the Hamilton- 
Jacobi differential equation in Kerr metric 

2.1 Determination of first geodesic integrals using the Hamilton- 
Jacobi approach assuming vanishing cosmological con- 
stant 

In this section we are going to review the method of integration of the Hamilton- 
Jacobi differential equation, when separation of variables is taking place and 
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apply it in deriving the timelike geodesic equations that describe the motion 
of a test particle in a Kerr spacetime. The more general case of Kerr-de Sit- 
ter spacetime is treated in section 6. The method for deriving the geodesies, 
assuming vanishing cosmological constant in the Kerr spacetime solution, was 
first applied from Carter [23]. 

The Hamilton- Jacobi differential equation is given by: 

dW / dW dW dW\ n 

-^ r +i?^;g 1 ,g 2 ,.-. ,g fc ;^-,^-,..- =0 (4) 

where if, the characteristic function, depends on the degrees of freedom <7i, <?2, • • • ,qk,t 
and the momenta pi := i = 1, • • • k. Also H, is a function of second degree in 
the partial derivatives of W. Assuming, that an integral W(t; qi, q<i, •• • , ol\, ■ ■ ■ ,a k ) 
has been found, where a\,- ■■ ,a k are integration constants, then the general 
integral of the differential equations of mechanics 

dqi dH 

dt dpi ' 

dpi dH 

dt dqi 

(5) 

is provided by 



— -P p - dW 

dai 1 ' dqi 

dW a dW 

da k dq k 



(6) 



where /?!,••• , /3 k denote new arbitrary constants. 

An interesting case, which was first studied in great detail by Stackel [25] 
and Levi-Civita [26] , is when one can integrate the Hamilton- Jacobi differential 
equation through separation of variables. In this case, the characteristic function 
W takes the form: 

n 

W = Y,Wi{qi) (7) 
i=i 

where every term on the right hand side of eq.(7) Wi(qi) depends only on one 
of the variables 

Assuming now, that the characteristic function separates in the Kerr metric 

W = W t {t) + W^) + W e + W r (8) 

where W$ , W r are functions of 9 and r coordinates respectively, and then plug 
this expression in the relativistic Hamiltonian in BL coordinates 

H(x a ,pp) = ±g' a 'p llPv (9) 
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and the relativistic form of eq.(4), S^^f^-f^ + M 2 — 0, one gets 

VK = -£ct + y ^dr + J V@d6 + L<j> (10) 



where 



and 



6 :=Q 



a=V - E ) + 



-A- cos 2 fl (ID 
sin 



i? := [(r 2 + a 2 )E — aL~\ 2 — A [^V + (L - s£) 2 + Q] (12) 

with A : = r 2 + a 2 — 2< ^y r . Also E, L are constants of integration associated with 
the isometries of the Kerr metric. Carter's constant of integration is denoted by 
Q. Differentiation of (10), with respect to the integration constants E,L,Q,[i 
leads to the following set of first-order equations of motion [23] : 

2 , 



2 cdt 

p lx 


r 2 + a 2 
A 


2 dr 
P dX 


= ±Vr 


2 de 

P dX 


= ±Ve 


2 d<p 
P dX 


A 



a(aE sin 2 6 - L) 



sin 



^ (13) 



where 



P:=E(r 2 +a 2 )-aL (14) 

We proceed now to discuss the exact solution of timelike geodesic equations 
(13) in closed analytic form, and apply them to the precise determination of 
Lense-Thirring precession of a satellite's spherical polar orbit in the gravitational 
field of Earth, as well as in stellar orbits aroung the galactic centre. We should 
mention at this point, the extreme black hole solutions a = 1 of spherical non- 
polar geodesies obtained in [24] expressed in terms of formal integrals. The more 
general case in the presence of the cosmological constant, will be discussed in 
section 6. 



3 Spherical Polar Geodesies 

Depending on whether or not the coordinate radius r is constant along a given 
timelike geodesic, the corresponding particle orbit is characterized as spherical 
or non-spherical respectively. In this section, we will concentrate on spherical 
polar orbits. 

Assuming zero cosmological constant, r = rf, where rj is a constant value 
and using the last two equations of (13) we obtain: 

d<j> _ 2£-aE + L/sm 2 9 

de ( ] 



G 



Equation (15) can be rewritten 

Maine ( a[ y+_iL a , L1 - aE + L/sin 2 e\ 

d<j> = 1 r " J (16) 

y/Q(l - cos 2 6) - [a 2 (^ 2 - E 2 ){1 - cos 2 6) + L 2 } cos 2 6» 

Now by denning z := cos 2 9, the previous equation can be written as follows: 

1 dz (aP L } 

dj> = —-—= x\ — -aE+ } (17) 

2 ^/z 3 a- z 2 (a + fl) + Qz I A 1-zJ 

where 

a:=aV-£ 2 ), f3:=Q + L 2 (18) 

It has been shown [27] that a necessary condition for an orbit to be polar 
(meaning to intersect the symmetry axis of the Kerr gravitational field) is the 
vanishing of the parameter L, i.e . L = 0. Assuming L = 0, in equation (17), we 
can transform it into the Weierstrafi form of an elliptic curve by the following 
substitution 

"—^sft (19 » 

Thus, we obtain the integral equation 

and this orbit integral can be inverted by the Weierstrafi modular Jacobi form 

3 

Z = p{<t>/A) (21) 

where A := -\ - aE^j , P> = E(r 2 + a 2 ). 
The Weierstrafi invariants are given by 

a 2 B 2 a[3 Qa 1 2 Qa 

92 = 12 + 12+T- — =12 (a + /3) -X 

a 3 a 2 /? c^Q _ _L_ 3 QoP_ Qa8 

93 ~ 216 + ~72~ + ~72~ + 216 ~ ~48 48~ ~ 216 ^ + ^ ' ~ ~48 48~ 

(22) 

We note that the discriminant A c := g\ — 27g 2 of the elliptic curve, vanishes 
for a = 0, 4 and this is consistent with the fact that circular orbits in the 
Schwarzschild spacetime do not have modular properties [3]. 
Now 

- E \ (23) 

3 For more information on the properties of the Weierstrafi function, the reader is referred 
to the monographs [34, 40], and the appendix of [10]. 
4 in this case a = 0, fi = Q + L 2 (= Q for polar orbits). 
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Then 



where 



# - dZ , A' (24) 

y/z*a! -z 2 (a' + 0) + Q'z 



* = r , g( ^l,-gUf-i 

I [£ + 1 



-EGMr 



i -i IGMr 
\ ± — — 2 — IT" 



a' 2 




a'/?' 


12 


+ Ty + 


6 


a' 3 
216 


a' 2 /?' 
+ 72 


+ 72 



4 

2 «' 3 „,/2n' n ,ia'rii 



(25) 



where a' := %,0 := Q/a 2 = Q'. Using Eq.(19) with a -» a',/3 -» /?', and 
integrating Eq.(24) we get 

c = p(0M') (26) 

with the Weierstrafi invariants > 53 given by the expressions 



g' = 1 -I £— + £ ^ (27) 

ya 216 72 72 216 48 48 y ' 

Equation (26) represents, the first exact solution of a spherical polar orbit as- 
suming zero cosmological constant, in closed analytic form, in terms of the 
Weierstrafi Jacobi modular form. 

Equation (26) can be rewritten in terms of the original variables as 

p(4> + e) = ^-cos 2 9-±(a" + 0') (28) 

where a" := a' /A 12 , 0' := /A 12 , Q" := Q'/A 12 with corresponding Weierstrafi 
invariants g'^g'l- Also e is a constant of integration. 



3.1 Frame dragging in spherical polar geodesies assuming 
vanishing cosmological constant 

We are now going to calculate the Lense-Thirring effect which as we discussed 
in the introduction is the imprint of rotating black holes or rotating mass distri- 
butions, i.e. the dragging of inertial frames. For this we shall need to calculate 
the periods of the Weierstrafi modular form. 

In order to calculate the periods of the Weierstrafi function we need to de- 
termine the roots of the cubic. The Weierstrafi invariants are given in terms of 
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the Kerr parameter a and the initial conditions E, Q as follows 5 : 



1 (a 2 (l-E 2 ) + Q) 2 1Q(1-E 2 ) 



12 a 4 A' 4 4 a 2 A' 4 

9 * = 432^[2« 6 (l-^ 2 ) 3 -3a 4 (l-i? 2 ) 2 g-3 a 2 (l-i? 2 )g 2 + 2Q 3 ] 

(29) 

while the discriminant A c of the cubic equation is given by the expression 

AC _ (i-£ 2 ) 2 Q 2 (« 2 (-l + £ 2 ) + Q) 2 .... 

256a* A' 12 (M) 

The sign of the discriminant A c determines the roots of the elliptic curve: A c > 
0, corresponds to three real roots while for A c < two roots are complex 
conjugates and the third is real. In the degenerate case A c = 0, (where at 
least two roots coincide) the elliptic curve becomes singular and the solution is 
not given by modular functions. The analytic expressions for the three roots of 
the cubic, which can be obtained by applying the algorithm of Tartaglia and 
Cardano [54], are given by 



ei = 



('2 



(q 2 (-l + ff 2 ) + 2Q) A 2 
12a 2 E 2 r 2 (GM/c 2 ) 2 
(2a 2 (-l + E 2 ) + Q) A 2 
12a 2 E 2 r 2 (GM/c 2 ) 2 
(a 2 (-l + E 2 )-Q) A 2 
12a 2 £; 2 r 2 (GM/c 2 ) 2 

(31) 



Since we are assuming spherical orbits, there are two conditions from the 
vanishing of the quartic polynomial R and its first derivative. Implementing, 
these two conditions expressions for the parameter E and Carter's constant Q 
are obtained, 

E 2 = 



{r 2 +a 2 )Z 

r(r 3 GM/c 2 + a 2 r 2 -3a 2 ^r + a 4 ) „ , 
Q - -i -2- J --a 2 E 2 

(32) 

where Z := r 3 + a 2 r - 3^fr 2 +a 2 ^f [27]. 

5 We have set jjt = 1. 
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Using the first and fourth line of Eq.(13) with L = 0, we obtain 



cdt r^2_p' - a 2 Esm 2 6 

~d4 ~ fP' - aE 

a 2 E cos 2 9 

= a + —pr^ h —ft, (33) 

Z^-aE ^-aE V ' 



or 

ct + £ = ad) + 



r 2 P'/A 4a 2 E/a" 



{-2aA'Y {-2aA') 

where we used Eq.(28) and the fact that, J p((f>)d(f> = —((4>), where ((z) denotes 
the WeierstraB zeta function. Also £ denotes a constant of integration. 

Our free parameters are a and r. For a given choice of our free parameters, 
E and Carter's constant Q are fixed with the aid of (32). 

The two half-periods w and uj' are given by the following Abelian integrals 
(for A c > 0) [53]: 

f°° dt , f £3 dt 

uj= / =, J = i\ (35) 

J ei - g 2 t - g 3 J-oo V-4t 3 + 32* + 53 

The values of the WeierstraB function at the half-periods 6 are the three roots of 
the cubic. For positive discriminant A c one half-period is real while the second 
is imaginary 7 . The period ratio is defined as r = — . 

After a complete oscillation in latitude, the angle of longitude, which de- 
termines the amount of dragging for the spherical polar orbit in the General 
Theory of Relativity (GTR) , increases by 

A</> GTR = 4lu (36) 

3.1.1 Dragging of satellite polar orbit in the gravitational field of 
rotating Earth 

In this subsection, assuming the Earth's gravitational field is described by a Kerr 
space-time, we apply the exact solution of the Kerr spherical polar geodesies 
obtained in the previous section, in determining the effect of Earth's rotation 
on the motion of its orbiting satellites. Thus we perform a precise calculation 
of the corresponding Lense-Thirring effect. 

Experimentally, a polar orbital geometry has the advantage that classical 
nodal precessions due to the even zonal harmonics of the gcopotential, are van- 
ishing. We note, that the Gravity Probe B (GP-B) mission, launched in April 

6 An alternative expression for the real half-period ui of the WeierstraB function is: ui = 

Vei - e3 f f1 (i §>b ffeff ), where F{a,(5, 1 ,x) is the hypergeomctric function 1 + + 

Q (q + l)/3(/3+l) 2 , 
1.2.7(7+1) 

We organize the roots as: ei > e2 > e$. 
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parameters 

= 329.432cm, r = r GP _ B = 7027.4 km 
329.432cm, r = r PO LARES = 8378 km 



half-period 

lu = 3.699746 x 1CT 11 
uj = 2.8421923 x 1CT 11 



predicted dragging 

A0 GTR = 0.164^ 
A</> GTR = 0.0969^ 



Table 1: Predictions for frame dragging of a satellite's polar spherical orbit in 
the gravitational field of Earth. The period ratio is r = 6.18642i. 



2004, retains a stricly polar orbital configuration during the science phase [28], 

8 9 

The Kerr parameter that corresponds to the angular momentum of the Earth 
is equal to a = 329.432 cm = 371.398(2GM0/c 2 ). For the radius of our 
spherical orbit we use the semi- major axis of the GP-B mission r = ?"gp-b = 
7027.4 km 10 . Then Carter's constant Q and the parameter E are determined 
by eq.(32). Then, the half-period u> has the value u> = 3.699746 x 10" 11 , which 
leads to the precise frame-dragging effect of A GTR = 0.164 arcsyr -1 . We 
repeated the analysis, for slight different values of radius r, and for fixed Kerr 
parameter. In particular we used as a value for r the semi-major axis of the 
Polares satellite, i.e. r = rp i a ros = 8378 km [31]. We list our results in table 1. 
We also compare our precise results with the Lense-Thirring (LT) formula [14]: 

A ^ LT - c2r 3 (1 2GJ e2)3 / 2 ( 37 ) 

in which J is the magnitude of the central body's angular momentum, r e the 
semi-major axis of the orbit and e the eccentricity of the orbit. For zero eccen- 
tricity, formula (37) predicts 0.164 arcs yr _1 and 0.096 arcs yr _1 for the radii of 
GP-B mission and Polares satellites, respectively. 

8 http:/ /cinstein. stanford.edu/ 

^Indeed if the orbit is exactly polar the precession due to Earth's quadrupolc moment 
J2, vanishes. Practically, it had originally been suggested by Van Patten and Everitt a two 
satellite experiment, with the two satellites in similar nearly polar orbits. The quadrupolc 
moment of the Earth, then produces precessions of opposite directions for the two satellites 
while the Lense-Thirring precession produces the same effect for both. To separate these two 
effects, the small angle between their orbits needs to be measured in addition to the precession 
rates of the two satellites. The distance of the two satellites at each polar crossing is measured 
by Doppler ranging. Substracting, the contribution of the quadrupole moment from the actual 
measurement determines the Lense-Thirring effect. Ciufolini [29] alternatively suggested, to 
make use of a satellite already in orbit LAGEOS and require only one new satellite LAGEOS 
II. Although the orbirt of LAGEOS, is far from polar, if the orbit of LAGEOS II is inclined 
at an exactly opposite angle with respect to the axis of the Earth, the Lense-Thirring effect 
can directly be measured, since the contributions from the quadrupolc moment cancel exactly. 
See also [33] for a different approach. 

10 We note the small eccentricity (0.014) of the GP-B satellite. 
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parameters 




half-period 




predicted dragging 


acaiactic = 0.52, r = 


III 


cj = 0.0515693 


A0 GTR = 


11.8188° per revolution = 42531.4 a !^- 

1 revolution 


acaiactic = 0.52, r 


50 


w = 0.00462023 


A^ GTR 


= 1.0589° per revolution = 3810 " c "- 

i revolution 



Table 2: Predictions for frame dragging from galactic black hole, with Kerr 
parameter acaiactic = 0.52 GM J^ H , for two different values of orbital radius. The 
values of the radii are in units of GMbh/c 2 The period ratios, t, are 2.92558i 
and 3.81583i respectively. 



parameters 

acaiactic = 0.9939, r = 10 
acaiactic = 0.9939, r = 50 



half-period 

ui = 0.0981121 
to = 0.00882902 



predicted dragging 



Ac/)' 
A<f> 



GTR 



GTR 



22.485° per revolution 
= 2.023°per revolution = 



= 80917.2 ! ' 
7281.66 



revolution 

arcs 



rnvolntion 



Table 3: Predictions for frame dragging from galactic black hole, with Kerr 
parameter acaiactic = 0.9939 GJ ^f ff , for two different values of orbital radius. 
The values of the radii are in units of GMbh/c 2 . The period ratios, r, are 2.5086 
i, 3.40336 i respectively. 



3.1.2 Precise calculation of dragging of stellar polar orbits from a 
galactic black hole 

Assuming, that the centre of the Milky Way is a black hole and that the structure 
of spacetime near the region Sgr A* , is described by the Kerr geometry as is 
indicated by observations Eq.(l), we determined the precise frame dragging 
(Lense-Thirring effect) of a stellar orbit with a spherical polar geometry. The 
results are displayed in table 2. 

We repeated the analysis for a value of the Kerr parameter as high as 
aGaiactic = 0.9939. Such high values for the angular momentum of the black 
hole, have been recently reported from X-ray flare analysis of the galactic cen- 
tre [30]. The results are given in table 3. 



4 Spherical geodesies with L ^ 

We now integrate Eq.(17) including the contribution from the parameter L. Let 
us define: 

n := / - ^ (38) 

J V4?-92t-gs 

thus £ = p(n+e), and A" f -. M = A" f —. & = A"U, A" := 

— \ — aE) and e is a constant of integration. Now 

1 T f dz 

~—L / (39) 

2 J (l-z)^z 3 a-z 2 (a + l3)+Qz 
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under the substitution (19) becomes 
La 



La 



53 



La f 
La r 



P'(n)dn 



(w - P (n + c )) v/4 P 3(n) - 32 P (n) - ff3 
rfn 

(«;-p(n + e)) 



La 



Log 



cr(n 



<r(n + e + w ) 



+ 2n c(v ) 



1 



(40) 

where w := j (l — = p(^o)- Also a(z) denotes the WeierstraB sigma 

function. Thus the equation for the orbit is given by 



= J d<j> = A"U - 



La 



Log 



cr(n 



«o) 



<j(n + e + v ) 



+ 2ncM 



(41) 



and p 2 '(v ) = 4p 3 (u ) - 32p(«o) - 93 = 4w 3 - g 2 w - g 3 . In terms of the 
integration constants, w is given by the expression: 



w = 



a 2 {l-E 2 ) a 2 (l- E 2 ) + Q + L 2 



12 



(42) 



Similarly using the first and third line of Eq.(13), we obtain for t the expres- 



sion 



r 2 + a 2 n an\ „ T , a 2 E^( la + /3 



+ ^>> (43, 



5 General solution for the time-like geodesies in 
Kerr metric 

In the general case, of non-spherical orbits one has to solve the differential 
equations 



dO 



r dr 

Vr 



where R(r) is a quartic polynomial given by 

R = [(r 2 + a 2 )E - aL~\ 2 — A [r 2 + (L - aE) 2 + Q] 



(44) 



(45) 



and 0(9) is given by equation (11). Note that this equation is an equation 
between two elliptic integrals, in the general case when all roots are distinct. 
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The left hand side can be transformed into an elliptic integral with variable z 
or £ in Weierstrafi normal form, see Eq.(19). In order to solve this differential 
equation and determine rasa function of 9, we will employ the theory of Abel 
for the transformation of elliptic functions [17]. 

The general case of non-spherical orbits is of importance for determining the 
precession of perihelion, periapsis, and perinigricon (point of closest approach 
to the black hole). 



5.1 Exact solution of the general Kerr geodesic using the 
transformation theory of elliptic functions 

Abel in [17] deals with the following question: 

" Find all possible cases in which the differential equation 

dy =±C , dX (46) 

^/(l-c 2 y 2 )(l-e 2 y 2 ) v /(l-c%2)(i_ e 2 a; 2 ) ' 

is satisfied by putting y an algebraic function of x, rational or irrational. " 

He explains that the problem may be reduced to the case in which y is a 
rational function of x. Below we shall describe his approach, which uses the 
idea of inversion of elliptic integrals and then we shall apply it to the solution 
of eq.(44). Abel's notation is as follows: From the elliptic integral 

* = / ^ (47) 

Jo ^{l-c 2 x 2 )(l-e 2 x 2 ) 

he defines the inverse function 

x := Atf (48) 

while 

w f 1/c dx J f 1/e dx 



r 1/c dx <j__ r 

2 7o v /(l-c 2 x 2 )(l-e 2 x 2 )' ~2 ~ J Q ^(l-c 2 x 2 )(l 



(49) 



and Afl = y(T— c 2 x 2 )(l — e 2 x 2 ). Also the solution of the equation Ai9' = Ai9 
is = (— l)" l + m %j + moj + m'ui' where m, m' are arbitrary integral numbers. 

Let y = ip{x) be the rational functions we are looking for and x = \i9,Xi — 
X'&i two roots of the equation y = ip(x), (y = ip(x) = ip(xi)). Let the radical 
on the left hand side of Eq.(46), be denoted by v^R, then 

-$= - ±Cd$ (50) 

Vr 

Then by changing x with xi, we have ±^L = ±Cdi?i, which gives d/di = ±d-d 
or after integration &i — e ± $ where e is a constant. The quantities \$ and 
A(i? + e) arc roots, thus we have 

y = ^(\$)=i>(\(v' + e)) (51) 
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Now 



y = $ (At?) = i> (A(t9 + e)) = V (A(t9 + 2c)) = • • • = V (A(0 + fee)) (52) 

where k denotes an arbitrary integral number. As the equation y = ip(x) has 
only a finite number of roots, there exist k and k' distinct such that A(t9 + fce) = 
A($ + k'e). Assuming n := k — k' > 0, and replacing § by d — fc'e, we get 

A(t9 + ne) = Ai? (53) 

Thus we get ■& + ne = (— \y n + m ' $ + TOW + m'u/, and since $ is a variable we get 
(-1)™+™' = 1 and 

m to' , . 

e = — w H J (54) 

n n 

where n :— mjn^n' := m'/n € Q, to + to' is even number. Assuming the degree 
of the equation y = ip(x) is such that it has roots other then A(t9 + fee), anyone 
of them has the form A($ + ei) where ei = /l^w + fjf^u)' , (/i l7 ^ G Q) and all the 
A(-i? + fee + fciei) are roots of the equation. Continuing in this fashion the roots 
of y = ip{x) are of the form 

x = A(i? + fciei + fc 2 e 2 + k 3 e 3 H h fc„e„) (55) 

where fci, &2, ■ ■ ■ fc„ are integers and the quantities ei, £2, • • • e„ are of the form 
+ fjfu)'. By denoting the values of the expression A(i9 + fciei + fc2£ 2 + ■ ■ ■ fc^Ci/) 
by At?, A(t? + ei), A(# + e 2 ), ■ ■ ■ A(i? + e m _i), and writing ^(2;) = |, where p, g are 
polynomials of degree to in ir, with no common divisor, the equation y — ip(x) 
is written 

p-qy=(.f- 9V){x - M)(x - A(0 + ei)) ■ • • (x - \(i9 + e ro _i)) (56) 

where the constants /, g are the dominant coefficients of p, q respectively. If /' 
and g' denote the coefficients of x m ~ x we have that 

/' - g'y = -(/ - gy) [A0 + A(0 + ei) + A(tf + e 2 ) + • • • A(t9 + e ro _!)] (57) 
and if we define 

ipd := Atf + A(t9 + ei) + X{& + e 2 ) + • • • A(0 + e m _i) (58) 
we get for y the expression 

V = -77 5" ( 59 ) 

By hypothesis, y is a rational function in x and so <p must be the same. The 
goal is to express ip rationally in x and to determine /, /', g, g' , ei, c\ and C in 
order that (46) is satisfied. 

First, we search for values of e which satisfy A(i9 — e) = \{d + e), that is to 
say A(i? + 2e) = Ai?. According to (54), 

6 = y w + t (60) 
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where m + m' e 2Z. In this case A(# + e) takes the following distinct values [17] 

At? = x, 
A(t? + w) = -Xd = -x, 



2 2 ec x 



1 



1 1 



3w u/ 1 



(61) 



All the other roots A(t? + e) come with a term A(t9 — e) and therefore we can 
write an expression for <ptf as follows [17]: 

iP & = \»+ k\(d + + k'\(d + | + y ) + fc"A(tf + ^ + y ) 

+A(0 + ei) + A(t9 - ei) + A(0 + e 2 ) + A(i? - e 2 ) + • • • + A(0 + e„) + A(t9 - e„) 



(62) 



where fc, fc', fe" are equal to or 1. 



5.1.1 Rational solution for y = ip(x), when k = k! = k" = 

In this particular case, we have 

ipd = Atf + A(tf + ei) + A(i?-ei) + A(i? + e 2 ) + A(i?-e 2 ) 

+ ■■■ + \(d + e n ) + \(d-e n ) (63) 

or written in terms of x 

Ae, 



tpd = x + 2x % —r 

4^ 1 - e 2 c 2 X 2 e,x 2 



(64) 



Thus, the prime condition that y is rational in x is fulfilled. Let S, 5' , 77, 77' be the 
values of 1? that correspond to y = 1/ci, — 1/ci, 1/ei, — 1/ei. Thus, for 3/ = 1/ci, 

1 — Ci w = = s'-^f'+jg-d^^pS 



or 



1 - ciy 



.9' - ci/' 



</?<5 



(65) 



where £ := g' + g<fd. Similarly, we get [17] 

g' + cif 



1 + ciy 
1 - eiy 
1 + eiy 



s 



tpil. 



(66) 
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From the expression of one gets 

x _ft = 1 + A lX + A 2 x 2 + ■ ■ ■ A 2n+1 x 2n+1 

<p5 (l-e 2 c 2 \ 2 e lX 2 )(l~e 2 c 2 \ 2 e 2 x 2 )---(l-e 2 c 2 \ 2 e n x 2 ) 1 ' 

For § = S the right hand side of (67) must vanish and similarly for d±e\ , ■ ■ ■ 5±e n 
with arbitrary 5. Thus we have that 



1 + A.. + ... + W" ^ (i-^) (i-ssttt 1 



1 - 



A(5 + ei)7V A(5-ei) 

X \ / X 



X(S + e n )J\ X(S-e n ) 



(68) 

Eq.(46) can be written as 

^(1 - c 2 y 2 )(l - e 2 ?/ 2 ) - ^ V(l - ^ 2 )(1 - e 2 ^ 2 ) (69) 

and we can see that the four functions l±cij/,l±ei?/ must vanish when x = 
±i,±i that is when = ±|,±^. Thus for 5 = |, J' = -|, jy = ^,7/ = 

we get 5 ' = Cl ./V (|) = ei/ V (Y) and 11 /' = £p (f ) = £p . A solution 
of this system is 

-'-"•F -^-^■-4) (70) 



and k is arbritrary. Then 



?y = (71) 



and 



12 



<Kf) ^ I A(f-ei)J I A(|-e 2 )J t A(f - e„) 

(72) 

where/, := (1 — e 2 c 2 \ 2 e 1 x 2 )(l — e 2 c 2 A 2 e 2 a; 2 ) • • • (1 - e 2 c 2 A 2 e„x 2 ). Similar expres- 
sions are obtained for 1 + , 1 ± -^r- . Then 1 - c 2 y 2 = (1 - ^|)(1 + 2&) = 

(1 _ c 2 x 2 )£, 1 - e 2 y 2 = (1 - e 2 x 2 )^ or 

yj(l - c 2 y 2 )(l - e 2 y 2 ) = ±^-y/(l - c 2 x 2 )(l - e 2 x 2 ) (73) 
n We have used the fact that Abel's function satisfies: (p (-tt) = — ¥>(^) >¥> 

12 Af = l/c. 
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where d := ^1 — x i (^-e 1 ) ) ' " {} ~ A^-e ) ) an< ^ a similar expression for d! by 



2 dy 



substituting w/2 — ► u//2. Then it is shown [17] that -ffi is a constant C and 
Eq.(46) is satisfied. 

Two other expressions for y are also provided [17] 



y = c 



x 1 - 



(1 - e 2 c 2 A 2 eiZ 2 )(l - e 2 c 2 A 2 e 2 x 2 ) ■••(!- e 2 c 2 A 2 e„x 2 ) 



and 



V=- (ec) 2n b Atf A( Cl + 0)A(ei -#)••• A(e„ + t?)A(e„ - 0) 



(74) 
(75) 



with 6 := A 2 eiA 2 e 2 A 2 e 3 • • • A 2 e„. For i? = w/2,i? = u//2, with corresponding 
values for y = 1/ci and 1/ei we obtain 

6 



1 

Cl 

1 

ei 



(_l)»^ e 2 « c 2 ™-i 



A^-eOA^-e,).-^ 



2n-l c 2n 



(76) 

As is discussed in [17] (see also [18]) the transformation found by Jacobi in [20], 
corresponds to ei = 2 2 " 1 , c = Ci = 1, and then e 2 = 2 n+i = -^l^n = ne\ and 



^A^ + ^A^-^-.-A^+^A^-i?) 



2n+l 2 



S A 



2n+l 2 



C = 



ei 



A (ot) a (jri) ' ' ' A pgr) 

t, / 1 u\ \ / 3 u>\ \ ( 2n-l u>\ 
. ^ 2n+l 2 y A \ 2n+l 2 J ' " A ^ 2n+l 2 J 



2n — 1 cj 
2n+l 2 

2 



V(l-y 2 )(l-efj/2) 



a 2n+l 



±C 



1 



2n + 1 2 



^(l-x 2 )(l-e 2 x 2 ) 



2n + 1 2 
= ±Cdtf 



2n — 1 uj 
2n+l2 



In addition in [17], two more cases are discussed i) those in which k = and 
k! or k" is equal to 1 and ii) k = 1. The general theorem for the first real 
transformation of arbitrary order n is also stated [17]. 



In Jacobi's language the indefinite elliptic integral x = J ■ 



V(i-V)(i-fcV) 

defines the elliptic function y = sn(x, fc).Thus, in Eq.(77), x — Xi3 — sn(?9, e) 13 . 

13 Onc can also switch to Lcgcndrc's notation by defining y = sin?/), x = simp'. Then the 



last line of (77) can be written as , a = 

J sin 2 i/> 

becomes negligible and then ip = sin -1 y = C J 



C 



d<p 



dip 



For large n the modulus ei 



(77) 
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5.2 Transforming the geodesic elliptic integrals into Abel's 
form 

The transformation 

* - es + (78) 
x z 

transforms the elliptic integral in Weierstrafi form into Abel's and Jacobi's form: 

f dx ■ 1 / dx (79) 

J ^A{x-e 1 ){x-e 2 ){x-ez) - e 3 J ^(1 - x 2 )(l - k 2 x 2 ) 

with k 2 = 

ei-e 3 

Similarly, the quartic can be brought to Jacobi's form y 2 = (l — x 2 ){l — k\x 2 ). 
Indeed, as we saw in section 3 we have 

dO 1 dz 1 d£ 

Ve 2^z*a-z 2 (a + f3)+Qz 2 ^ _ g ^ _ g3 

where 

92 = ±(a 2 (l-E 2 ) + Q + L 2 ) 2 -^(a 2 (l-E 2 )) 

93 = ^(a\l-E 2 ) + Q + L 2 ) 3 -^(l-E 2 ) 2 -^(a 2 (l-E 2 )(Q + L 2 )) 

(81) 

In this case the three roots of the cubic 4£ 3 — g 2 ^~ 93 are given by the expressions 
ei = ^ (a 2 (l - E 2 ) + L 2 + Q + 3^/(a 2 (-l + E 2 ) + Q) 2 + 2a 2 L 2 - 2a 2 E 2 L 2 + L 4 + 2L 2 q) 
e 2 = ^ (a 2 (l - £ 2 ) + L 2 + Q - 3^(a 2 {-l + E 2 ) + Q) 2 + 2a 2 L 2 - 2a 2 E 2 L 2 + L 4 + 2L 2 Q) 
e 3 = l(a 2 (E 2 -l)-L 2 -Q) 



(82) 



Thus the Jacobi modulus k 2 — 62 63 is given by 

fc2 = -a 2 (-l + E 2 )+L 2 + Q- y/a" 4 (-l + E 2 ) 2 - 2a 2 {-\ + E 2 )(£ 2 - Q) + (£ 2 + Q) 2 

-a 2 (-l + £ 2 ) + L 2 + Q + v/a 4 (-l + -E 2 ) 2 - 2a 2 (-l + E 2 )(L 2 - Q) + (L 2 + Q) 2 

(83) 

It has the correct limit for polar geodesies with L = 

k\L = *) = a ^^l (84) 

Also 

1 2V2 



Vei 63 J-a?(-l + E 2 ) + L 2 +Q + y/a*{-l + £ 2 ) 2 - 2a 2 (-l + E 2 ){L 2 - Q) + {L 2 T~Qf 

(85) 
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and it also has the correct limit for polar geodesies 

1 2 



Thus, we get 



(86) 



["=- 1 -- r J = [ - d l (87) 

7 VQ 2 Vei - e 3 J J( 1 _g2 ) ( 1 _ k2 g2 ) 

Now we shall discuss the transformation of the right hand side of Eq.(44). 
Inspired by the work of Jacobi and Abel, Luchterhand obtained the following 
transformation formula [21]: 

- dx = 9y (88) 

MVfi-^fi-^) V(y-a)(y-P)(y-i)(y-S) 

where the Jacobi modulus fci and the coefficient M, are given in terms of the 
roots a, (3, 7, S of the quartic, by the following expressions 



ki= y/ , i r S) !! i ~t M= y /{a- 1 ){fi-S) (89) 



v/(a- 7 )(/3-<5) 
and the integration variables are related by 

1 - x (7 - S)(y - a)(y - f3) 



(90) 



1 + x {a - 0){y - 7 )(y - J) 

Also we assume that the roots a, /3, 7, <5 of the quartic are real and are organized 
in the following ascending order of magnitude: a > [3 > 7 > 5. 

At this stage, we have succeeded in transforming both sides of Eq.( 44), into 
Abel's form. We can also provide a nice formula for r in terms of the Jacobi's 
sinus amplitudinus function 14 

( 7 -J)(r- a )(r-/?) = +-sn (m/ j§, fr) 
(a-/?)(r-7)(r-<y) 1 + S n (m/.*|, 

A systematic investigation of the phcnomenological aspects of the above 
exact solution like the periapsis precession will be reported elsewhere. 



6 Separability of Hamilton-Jacobi's differential 
equation in Kerr-(anti) de Sitter metric and 
derivation of geodesies. 

In the presence of the cosmological constant we proved (see Appendix A for de- 
tails) the important result that the Hamilton-Jacobi differential equation can be 

14 wc wright r f dr = [ — ; — 2s = r -41. 

& J VR J ^(r- a )(r-t3)(r-~/)(r~6) J M 7 (1 - x * ) (1 - h\ ^ ) J V& 
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solved by separation of variables. Thus in this case, the characteristic function 
separates and takes the form 



y/[Q +(L- aE) 2 - 2 jgg cos 2 fl] A 9 - gggg ^ 
W = -Set + L<j> + / <i# 



S 2 [{r 2 + a 2 )E - aL} 2 - A r (^ 2 r 2 + Q+ ~ 2 (L - aE) 2 ) 

dr 

By differentiating now with respect to constants of integration, Q, L, E, /i, 
we obtain the following set of geodesic differential equations 

dr f d6 



IB! J V& 

P 2 ^ = ~ A (aE sin 2 e - L) + ^ [(r 2 + a 2 )E - ail 

rfA A sin 2 (9 v ; A r LV ; J 

2 _ S 2 (r 2 + a 2 ) [(r 2 + a 2 )E - aL] a~ 2 (aEsm 2 - L) 

° P dX ~ A r A g 



t£ - (92) 



where 



2\ 



J?' := E 2 [(r 2 +a 2 )E- aL] 2 - A r (^ 2 r 2 + Q + E 2 (L- aE) 

& := [Q + (L - aE) 2 E 2 - ^ 2 a 2 cos 2 9} A e - E 2 (aE ^ ~ ^ (93) 

sin 6 

The first line of Eq.(92) is a differential equation that relates a hyperelliptic 
Abelian integral to an elliptic integral which is the generalisation of the theory 
of transformation of elliptic functions discussed in section 5, in the case of non- 
zero cosmological constant. The mathematical treatment of such a relationship 
was first discussed by Abel in [19]. 

6.1 Exact solution of spherical polar orbits with a cosmo- 
logical constant 

Using the second and the fifth line of Eq.(92), for L = and assuming a constant 
value for r, wc obtain 

HVE , a~?(r 2 +a 2 )E 



d4> ~^A^ + 

R A„ 

(94) 



-E 2 Ea + BAg 
A e V& 
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where B 



Similarly, using the third and fifth line we obtain 



cdt E 2 (r 2 +a 2 ) 2 E a 2 E 2 Esm 2 6 

(95) 



d6 ~~ A r V& AeVO 7 
rA -a 2 E 2 Esm 2 e 



A e VW 

— 2 t 2 , 2\2 rr 

and T := a (r g } £ . 

Now using the variable z = cos 2 9, we obtain the following system of integral 
equations: 



[ Z 2 aE/2 dz f B(l + 4M/(-2) 



•/ v/0) v/M 

or 

(«i + P 1 z)dz 

c t = rkh+t^L (97) 

where /(z) = z(l-z)(Q + z(Qa 2 A + S 3 a 2 S 2 - A i 2 a 2 ) + z 2 (- M 2 a 4 A))(l + a 2 | x )2. 
Also we have defined 



7i 



„ 2 a£ _ l aS 2 (r 2 + a 2 )g 

"2 2 " A r 

l a5 2 (r 2 + a 2 )£a 2 A 

2 A r 3 

lS 2 (r 2 +a 2 )£ a 2 S 2 E 

1 

2 A r 2 



a 2 AS 2 (r 2 +a 2 ) 2 £ a 2 E 2 E , s 

Sl = —6 a; — (98) 

Equation (96) is a system of equations of Abclian integrals, whose inversion 
in principle, involves genus-2 Abelian-Siegelsche modular functions. Indeed, this 
system is a particular case of Jacobi's inversion problem of hyperelliptic Abelian 
integrals of genus 2 [44]- [50] (see Appendix B for details). Then, one can express 
z as a single valued genus two Abelian theta function with argumenents t, </>. 
However, since the polynonial f(z) of sixth degree posses a double root it may 
well be that the Abelian genus-2 theta function degenerates and the final result 
can be expressed in terms of genus-1 modular functions. This issue will be 
investigated in detail elsewhere. 
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6.2 Equatorial geodesies including the contribution of the 
cosmological constant 

The equatorial geodesies (i.e. 9 — n/2,Q = 0), with a nonzero cosmological 
constant, may be obtained by Eq.(92) for the particular values of Q,8. The 
characteristic function in this case, has the form 15 

W = -Ect + J ^~dr + L<f> (99) 
and the geodesies are given by the expressions: 

dr dX 



2# a{l + \a 2 A) 2 (E{r 2 +a 2 )-La) , 1 2 2 

r dX = (l-A^+^,2^ +(L-aE)(l + -aA) 

2 dt (l + ±a 2 A) 2 (r 2 + a 2 )[(r 2 +a 2 )E-aL] l 2k s2 fT m 

cr 2 — = - 2 ^11 >. L + (1 + - a 2 A) 2 a(L - aE) 

(100) 



where 

R' = (1 + ia 2 A) 2 [((r 2 + a 2 )E - aL) 2 - A r ((L - aE) 2 )] (101) 
for null-geodesics and 

R' = (1 + \a 2 A) 2 [((r 2 + a 2 )E - aL) 2 - A r ((L - aE) 2 )] - A r (p 2 r 2 ) (102) 

for time-like geodesies. 

6.3 Various limits of equatorial geodesic equations. 

For A = 0,a^0 (Kerr limit) Eqs(100) reduces to 

dr dX 
VR ~ r 2 

2 dcj> a(E(r 2 + a 2 ) - La) 

r dX r 2 +a 2_2_GMr- + ( L - aE ) 

2 dt (r 2 +a 2 )[(r 2C +a 2 )E- a L] 

r c— = 1 + a{L — aE) 103 

dX A 



15 See appendix A for details. 
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Similarly for a = 0, A ^ (Schwarzschild-de Sitter limit), Eqs(lOO) reduce to 

dr dX 



J[r*E* - ((1 - ±£y - 2^H) [ M 2 r 2 + V]} 



r 2 



r Tx = L 

E 



3 rc 2 

(104) 



The exact solution of (104) have been given in [3]. 



6.4 Circular geodesies with a cosmological constant 

In this section we shall study the circular equatorial geodesies and we shall prove 
the result: 



dt _ rj 3 ± [(rj 3 ) 2 - r^rjo] 



1/2 



00 



a 1 -y/r 
c c 



/ / GM Ar \ 
V l^T 7 3 / 



(105) 



where r* fc denotes the Christoffel symbols. 

The way to arrive to this result is to demand Tj k u : 'u k = 0, together with 
the conditions: 

1 dr 

u 1 = — = 
ds 

2 d6 n 
u 2 = — = 

ds 

(106) 

Alternatively since we have derived the general form for equatorial geodesies 
by requiring m 1 = or equivalently that R'(r) — (and its first derivative also 
vanishes) will lead us to a condition for eliminating the constants of integration 
and we should arrive in (105). 



We also derive the following expression for u 3 = 4^. 



3 d<j) 
ds 



, dt 2 dt 
ff00( ^ ) + 2903 ^ + 933 . 



-1/2 

(107) 
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This is equivalent to the following equation 



ds 



^ (1 _A (r2 + a2) _ 2 GM ) (-)% 2 |^. 



ca(r 2j ra 2 ) -. 

(i+4^) 2 



dl 



(2GMr)/c 2 a 2 (r 2 +q 2 ) 
r 2 H 2 + , , a 2 A 

1_t ~ 3 



(108) 



[36]. 



The above expression may be obtained from the metric 

ds 2 = g QO dt 2 + 9ll dr 2 + g 22 d9 2 + g 33 d<jy 2 + 2g 03 d(f>dt (109) 
Some aspects of the equatorial circular geodesies have been investigated in 



7 Conclusions 

In this work, we have investigated the motion of a test particle in the gravita- 
tional field of Kerr space-time with and without the cosmological constant. The 
geodesic equations were derived by integrating the relativistic Hamilton-Jacobi 
partial differential equation by separation of variables. As we saw in the main 
body of the paper, the structure of the exact solutions of the corresponding 
geodesic equations is quite rich and mathematically very interesting. The in- 
tegration of the geodesic equations involves the inversion problem of Abel and 
Jacobi for hyperelliptic and elliptic Abclian integrals, as well as the theory of 
transformations of elliptic functions. As we saw, in the case of spherical polar 
orbits with a cosmological constant the solution involves directly the solution of 
Jacobi's inversion problem for a system of Abelian integrals. 

Physically, the new types of relativistic motions lead to the interesting phe- 
nomenon of frame dragging of inertial frames. In this respect, we have applied 
the solution of spherical polar geodesies in Kerr space-time in two situations. 
First, by modelling the gravitational field of Earth by the Kerr geometry, we 
determined the longitudinal dragging of a satellite's polar spherical orbit around 
the Earth using as radius the semi-major axis of the polar orbit of GP-B mission 
launched in April 2004. Secondly, by assuming the galactic centre is a rotating 
black hole, according to the interpretations of recent observations of the galactic 
centre [22], whose Kerr parameter is determined by experiment as for instance 
in Eq.(l), we determined the longitudinal dragging of a star in spherical polar 
orbit around and close to the galactic centre. In principle, these predictions can 
be tested by observations. 

A systematic investigation of the physical implications of the exact solutions 
of the geodesic equations for the various possible types of motion for the test 
particle, in the presence of the cosmological constant, will be reported elsewhere 
[52]. For instance, one can investigate the effect of the rotation of the Sun in 
the perihelion precession of Mercury and a comparison can be made with the 
cosmological constant contribution that has been calculated in [3]. 

It can also be applied in determining the cosmological constant effect in 
rotating galaxies, and in particular in the velocity rotation curves. Assuming 
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the galactic centre is a rotating black hole with the Kerr parameter measured as 
in eq.(l), what is the effect of the cosmological constant in the motion of stars 
as a function of the distance from the galactic centre ? 

Probing general relativity in the strong-field regime will certainly be one of 
the most exciting endeavours. 
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A Separability of Hamilton-Jacobi equation in 
the presence of a cosmological constant 

In this appendix, we will prove that the Hamilton-Jacobi diffential equation pre- 
serves the important property of separability in the presence of the Cosmological 
constant, which has been used in deriving the geodesic equations in section 6 
from the characteristic function W. 

The non-zero elements of the inverse metric are 



9°° 



9° 3 



9 11 



9 22 



S 2 [(r 2 + a 2 ) 2 E - a 2 sin 9 2 [(r 2 + a 2 )E - 

c 2 A r A eP 2 
a~ 2 (6GMr + c 2 A(a 2 +r 2 )p 2 ) 
3c 2 A r p 2 A e c 



P 2 
Ae_ 

P 2 



^ = S2 ^-^J (U0) 
The Hamilton-Jacobi differential equation takes the form: 

[{r 2 + a 2 ) 2 - - a 2 sin 2 9[{r 2 + a 2 )E - 2^£l]] E 2 sin 2 9 S 2 (sin 2 6A e a 2 - A r 



2 



p 2 A r A e sin 2 9 p 2 A r A g sm z 9 

2aZ 2 (m^ + A p 2 (a 2 + r 2))(_£L) sin 2 e 

p 2 A r A e sin 2 9 

A r fdW r \ 2 -A e (dW e \ 2 , 2 



+ 
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or 

— [(r 2 + a 2 )£ - ah] - a?{L - aE) 2 — + a?(L- aE) 2 — - — - — '- + 

A r A r Ag Ag sin a 

2 / ™„ \ 2 



+<'V + <-^) +<-»>(^) =» 

(112) 



For equatorial geodesies the Hamilton- Jacobi differential equation takes the 
form 

g 00 (-Ec) 2 + g^L 2 + ' + 2g™(-EcL) + /i 2 = (113) 

Using eq.(llO) with 9 = tt/2 in (113) we obtain 



dW r _ y/E 2 [{{r 2 + a 2 )E - aL) 2 - A r (L - aE) 2 } - A rf i 2 r 2 



dr 



(114) 



B Definitions of genus-2 theta functions that 
solve Jacobi's inversion problem 

Riemann's theta function [49] for genus g is defined as follows: 

9(u) := e 2wiun+iwQn2 (115) 



t»l,- , n g 



where fin 2 := flun 2 + • • • 2£lyiri\n2 + • • • and un := u\ni + ■ ■ ■ u g n g . The 
symmetric g x g complex matrix Q whose imaginary part is positive definite is 
a member of the set called Siegel upper-half-space denoted as £s g ■ It is clearly 
the generalization of the ratio of half-periods t in the genus g — 1 case. For 
genus g = 2 the Riemann theta function can be written in matrix form: 



0(u,fi) = e mtnnn+2n 

7T*( m ri2 ) 
e V 



£2 2 

Ul 



e 

111,112 



/ On 


0l2 \ 




1 fii 2 


«22 i 


(:) 



+27ril rii n 2 



«2 



(116) 

Riemann's theta function with characteristics is defined by: 

@(u;q, q') := e 27riti(ri+ 9 ')+i7rf2(n+<j') 2 +27ri<j(ri+(2') (117) 
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herein q denotes the set of g quantities qi, ■ ■ • ,q g and q' denotes the set of g 
quantities q[, ■ ■ ■ , q' g . Eq.(117) can be rewritten in a suggestive matrix form: 



e 



q 
q 



(u.fi) = e m ^ n+q '^ {n+q '^ +2m ^ n+q '^ u+q \ q.q'eQ 9 (118) 



neZs 



The theta functions whose quotients provide a solution to Abel-Jacobi's 
inversion problem are defined as follows [50]: 

9(u;q,q) := ^ e au 2 +2hu(n+q')+b(n+q') 2 +2i-Kq(n+q') (119) 

where the summation extends to all positive and negative integer values of 
the g integers n\, ■ ■ ■ ,n g , a is any symmetrical matrix whatever of g rows and 
columns, h is any matrix whatever of g rows and columns, in general not sym- 
metrical, b is any symmetrical matrix whatever of g rows and columns, such 
that the real part of the quadratic form bm 2 is necessarily negative for all real 
values of the quantities mi, • • • , m g , other than zero, and q,q' constitute the 
characteristics of the function. The matrix b depends on \g{g + 1) independent 
constants; if we put iiril = b and denote the (/-quantities hu by iirU, we obtain 
the relation with Ricmann's theta function: 

e(u;q,qi)=e au2 e(U;q,q') (120) 

The dependence of genus-2 theta functions on two complex variables is de- 
noted by: 6(u;q,q') = 9(ui,u 2 ;q,q'), the dependence on the Siegel moduli 
matrix f2 by: 9(ui,u 2 ,Q;q,q'). To every half-period one can associate a set of 

characteristics. For instance, the period u a - ai = | ^ j jj ^ while 9{u) is a 

theta function of two variables with zero characteristics, i.e. 9{u) — 9(u; 0, 0) = 

|j H (u, Ct). Also, Weierstrafi had associated a symbol for each of the six 

odd theta functions with characteristics and the ten even theta functions of 
genus two. For example, 9(u) is associated with the Weierstrafi symbol 5 or 
occasionaly the number appears as a subscript, i.e. 9{u)$. 

Let the genus g Riemann hypcrclliptic surface be described by the equation: 

y 2 = 4(x -ai)---(x- a g )(x - c)(x - a) ■ ■ ■ (x - c g ) (121) 

For g = 2 the above hyperelliptic Riemann algebraic equation reduces to: 

y 2 — 4(x — a\)(x — a 2 ){x — c)(x — c\){x — c 2 ) (122) 

where a\, a 2 , c Ci, c 2 denote the finite branch points of the surface. 

The Jacobi's inversion problem involves finding the solutions, for Xi in terms 
of Ui, for the following system of equations of Abelian integrals [50]: 

U x + ■ ■ ■ + Ui 9 =U\ 

: + ••• + : : 
u x g ^ + ■■■ + u x g ^ a ^=u g (123) 
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where u x,>1 



x dx ,x,n 

V ' v 



C x xdx 

Jfj, y ' ' 



x x 9 1 dx 

y 



i jil y ' * J fj, y > y J fj, y 

For 5 = 2 the above system of equations takes the form: 



f 



r xi dx 


r X2 dx 


J ai v 


Ja 2 y 


1 x dx 


f X2 x dx 






V 


/o2 V 
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= u 2 



(124) 



where Ui,U2 are arbitrary. The solution is given by the five equations [50] 
9 2 {u\u b ' a ) 



= A(b - xi)(b - x 2 ) ■ ■ ■ (b - x g ) 

= A(b - xi){b - x 2 ) 
(b-x\)(b - x 2 ) 



= ±- 



y/e* iPP 'f'(b) 



(125) 



where f(x) = (x — ai)(x — a 2 )(x — c)(x — ci)(x — c 2 ) 1 and e = ±1 according 
as u b ' a is an odd or even half-period. Also 6 denotes a finite branch point and 
the branch place a being at infinity [50]. The symbol Q(u\u h ' a ) denotes a genus 2 
theta function with characteristics: 0(u;q,q') [50], where u,— (ui,u 2 ), denotes 
two independent variables, see appendix A for further details. From any 2 of 
these equations, eq.(125), the upper integration bounds x\,x 2 of the system of 
differential equations eq.(124) can be expressed as single valued functions of the 
arbitrary arguments u\,u 2 . For instance, 



xi = Ol + 



\u\u ai > a ) 



and 



X2 



Ai(ar 2 -ai) 9 2 (u) 

(0 2 - Ol)(0 2 + Ol) + — % ffi(u) g2(M) 



(126) 



2(ai - a 2 ) 



± 



where 



(a 2 - ai)(o 2 + ai) + -fc- 



2 (lt|« a l- a ) 



8Hu) 



1 6> 2 ( t t|M°2.°) 

a 2 e 2 !^) 



4(ai - a 2 )rj 



2(ai - a 2 ) 



. a 2 6» 2 (du ai > a ) ai 6> 2 (du a2 ' a ) 
ry := a 2 ai(ai - a 2 ) ; :T — h 



Ai 6» 2 (m) A 2 6 2 (u) 



(127) 
(128) 



Also, A; = ±- 



Ve" ipp 7'(ai)' 

The solution can be reexpressed in terms of generalized Weierstrafi functions: 



(12) P2,2(u) ± yj Pl 2 {u) + 4p 2 ,i(w) 

c fc = o ! fc = 1, 



(129) 
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where 

P2,2(«) = " " (130) 

di — & 2 

and 

/ \ a, e 2 (tt|« a 2- a ) a , e 2 (tt|« a i.°) 

«,,(») - "*»"'■ : " 2) : + " dsn 

ai - a 2 

Thus, xi,a; 2 , that solve Jacobi's inversion problem Eq.(124), are solutions of a 
quadratic equation [45, 50] 

Ux 2 - U'x + U" = (132) 

where U, U', U" are functions of u\,U2- In the particular case that the coefficient 
of x 5 in the quintic polynomial is equal to 4, U = 1, U' = p 2 , 2 (u), U" = p 2 ,i(u)- 
The matrix elements fty , f2 y can be explicitly written in terms of the half- 
periods U?> a . For clarity, U^ £3 = Q xdx/y.U^ 3 = Q dx/y etc. The 
roots have been arranged in ascending order of magnitude and are denoted by 
e2 S , e2g-i, ■ ■ ■ ,e ,g = 2, so that e 2j ,e 2j _i are respectively c g - i+ i, a g - i+ i and 
eo is c. For instance, the matrix element hn = „ 4 ,<i 3 ^ ,1\ e2 , ei e4 , e3 , x 7ri, 



while On = 



2(t/r 

c/ ei,eo c/ e 4 ,e 3 _ c/ ei,eo !7 e4,e3 



Eq.(96) has the correct limit for vanishing cosmological constant. Indeed, 
for A = 0, P 1 = and 

(133) 

v /z(l-z)(Q + za 2 ( J B 2 -l)) 

which can be transformed to Jacobi's form by the substitution z = y 2 . Then 
the solution is given by the Weierstrafi function of section 3. 

C Transformation Theory of Elliptic functions 
and Modular equations 

One of the applications supplied by the transformation theory of Elliptic func- 
tions, which is of great importance in Number theory, are the modular equations 
described below [17, 18, 20]. 

For a rational solution of the differential equation 

d y = c dx fl 34) 

^{l-y 2 ){l-e 2 y*) ^(l-^Xl-e^) 

the necessary conditions among the periods 

K{e 1 ) = a CK(e)+a 1 CiK , (e) 
iK{ex) = b CK(e) + b 1 CiK'(e) 

(135) 
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with the period ratios (moduli) of the associated modular theta functions being 
given by 

r = b0 + hT ' (136) 
a + air' 

are also sufficient, when 

a bi — ai&o = n (137) 

is a positive integer number. The integer n is called the degree of transformation. 

Equation (137) for oo,&i,ai,6o <= Z when viewed as the determinant of a 
matrix £ GL(2, Z), sometimes is called a modular correspondence of level n. 

It can be shown that the inequivalent reduced forms of modular correspon- 
dences 

a oi 
bo h 

are of the form 





16£ q' 

where q a positive part of n represents, q' :— ^, and < £ < q' — 1. For instance 
for n = p a prime number, there are p + 1 inequivalent reduced forms of the 
form 16 

1 \ / 1 \ / 1 \ / 1 \ / p 
p ) \ 16 pJ'V 16.2 p ) ' " ' \ 16(p - 1) P J ' V 1 

Also the multiplication factor C in Eq.(134) is given by 

q K(e) 

which for a degree of transformation that is a prime number (n = p) is equal 
to K ^ll) or (l/ n ) • Eq.(138) can be rcexprcsscd in terms of Jacobi theta 
functions as follows 

C= l - (139) 

The modular equations are equations relating the Jacobi modulus e(pr) to 
e(r) which are of the form 




(140) 



16 



In a more familiar notation these classes of inequivalcnt reduced forms are 
q ^ J , ad = n, (a, b, d) = 1, < b < d and a,b,d g Z 
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where denotes the Legendre symbol 17 . Equivalently modular equations 

can be written in terms of the absolute modular invariant function j(r), and 
relate the reduced absolute modular invariant j* to j by polynomial equations 
of the form 

WJ) = (141) 

where j* :— j.cti = j ( aT £ b ) The explicit form of <I> 2 (j * , j) = 0, has been given 
in [41]. 

An additional example considered in [18] is that of the differential equation 

dy = ryrr . dx (142) 

v/(l-y 2 )(l + eV) v/(l-x 2 )(l + e 2 x 2 ) 

where y = ±ie n x jj^^^^y^^^^^T) . while thc modulus e is given 

u v 1 _ r n+l( ( uj \ _ _ _ / (2n-l)h) xx2 j C — { y("/( 2 "+ 1 ))-'V : '(»"/( 2 "+ 1 )) ) 1 

byl-e W 2(2n+1) J ^ 2 (2n+i)^ andG-^ - ) ... < ^ ( |^i f } J e - 

For n = 1 it gives e = \/3 + 2, <^ 2 (f ) = 2\/3 - 3. Also here ifd = x, where 
i 1 ) = f = 

J ^/(l-a; 2 )(l+e 2 a; 2 ) ' 

D Quartic and cubic for spherical orbits assum- 
ing A = 

In the case of spherical orbits (polar with L = 0, or non polar) and with the 
assumption of a vanishing cosmological constant, the test particle's radial co- 
ordinate will be stable at some value r/ if R(r) vanishes at r = rf and goes 
negative nearby. This will be the case if 



R(r f ) = 0, 

(143) 



— — r = Tt = 

dr 



and the second derivative of thc polynomial R(r) negative. 

These conditions result in the solution of a quartic and a cubic algebraic 
equation. Let us describe first how we can solve the quartic equation: 

For a quartic polynomial Ferrara (1545) solved the general equation of degree 

4 : 

Pi{x) = x 4 - ax 3 + c 2 x 2 - c 3 x + c 4 = (144) 

as follows: using the substitution x — > x + Ci/4, we can eliminate the cubic term 
(ci = 0). The reduced equation is equivalent to 

(x 2 + c ) 2 = yx 2 + c 3 x + Co - c 4 (145) 

17 (4)= e4 ^- 
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where cq = ^(c 2 + y). The discriminant A of the right hand side (RHS) is: 



A = 4- 4y(c 2 - c 4 ) = 4 - 4y (^(c 2 + yf - c 4 j (146) 

Its vanishing results in a cubic in y to be solved. Then the RHS has a double 
root x a = — and x 2 + c = ± v ^;(x — xq) maybe solved by the extraction of 
one more square root. 

As we described, the vanishing of the discriminant A = results in a cubic 
equation: 

-y 3 - 2c 2 y 2 + y(4c 4 - c 2 ) + c 2 = (147) 
or y 3 + 2c 2 y 2 + y{c\ — 4c 4 ) — c 2 = 0. Under the substitution 

v = -^- 2 -f ( 148 ) 



we get the cubic 



-64£ 3 + ^4 + 16c 4 ) - A c 3 _ c 2 + 8^4 (14g) 



or £ 3 + £m — n = with 



m , 4 2 
m := - — , m = -c 2 + 16c 4 , 



n ' i 2 , 2 8c 2 c 4 

" : =777 »« =-^ c 2-^H (150) 

64 27 2 3 3 

Then the roots are: 

For polar orbits the coefficients c 2 ,C3,c 4 in terms of the constants of inte- 
gration E, Q, (L = 0) and the Kerr parameter a are given by 

- 2 -^(Q + a 2 E 2 ) 

C3 = - 



c 2 = a 



E 2 - 1 

2 Q 

E 2 -l' 



c 4 = —gr— x ( 151 ) 
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